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Abstract 

The method recently proposed by Skala and Cizek for calculating perturba- 
tion energies in a strict sense is ambiguous because it is expressed as a ratio of 
two quantities which are separately divergent. Even though this ratio comes 
out finite and gives the correct perturbation energies, the calculational pro- 
cess must be regularized to be justified. We examine one possible method of 
regularization and show that the proposed method gives traditional quantum 
mechanics results. 
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Recently, in a letter in this journal |jT[, Skala and Cizek (SC) proposed a method to 
calculate perturbation energies using non-sqaure-integrable functions. The method of SC is 
further augmented in a comment by Guardiola and Ros (GR) p[. The purpose of our present 
comment is to further point out that in a strict sense, the SC method for the perturbation 
energies can result in a ratio of two divergent quantities, so that while this ratio may remain 
finite in a practical numerical calculation, a regularization procedure is needed to justify 
the finite result. We have examined one such possible regularization procedure and made 
contact between the SC method and traditional quantum mechanics (QM) results. 

Briefly, the SC method regards the n-th order perturbation equation as a parametric 
differential equation with En as the parameter, 

{Ho - Eo)MEn, x) = [En - V;)^o(a;), (1) 



where 



/^o = -|j + Vi, (2) 



and 



n-l 



K^O = Vilpn-l - J2 Ei^n~i- (3) 

i=l 

Here we have chosen to introduce Vn as the effective perturbation in the n-th order equation. 
For n = 1, the sum in R.H.S. of Eq. (|]) vanishes and Vi is the same as the real perturba- 
tion Vi. We also adopt the convention that all the wavefunctions are physical unless the 
dependence on the parameter En is explicitly displayed. For this form of Eq. (|l]), the n-th 
order equation is similar in form to the first order equation. V^ is a known function since all 
lower order quantities are assumed known. In traditional QM perturbation theory, one left 
multiplies Eq. (||) with the zeroth order wave function ipQ. Under the usual normalization 
conditions, 

(^o|^o) = l and (^o|^i) = Vz ^ 0, (4) 



one recovers the standard QM result 

En = (V-OIKI^O) = {MVl\^n-l). (5) 

Once En is correctly obtained, ipniEnjX) is obtained by solving the ordinary differential 
equation in Eq. (|I|). 

SC propose that instead of obtaining En in the standard way first, one treats Eq. (0) as 
a parametric ordinary differential equation with En a parameter, and go on to show that 

^ ^ -^n(0,a^o) ^ -■^n{0,Xo) ,g. 

lljn{l,Xo) -1pn{0,Xo) ~ F{Xo) 

where xq is a point such that the boundary conditions 

i)n{En,Xo) =0, n = 0, 1,... (7) 

are met for the physicsal energies En- Since neither nor 1 is the necessary n-th order energy 
correction En, the functions ipn{0,x) and ?/'„(!, x) are in general not square integrable and 
hence the name for the method. 

Since the ground state wave function vanishes only at the end points of the boundary^ 
and the nodal points of the wave functions of the excited states shift upon turning on the 
perturbation, the only choice for xq consistent with the boundary conditions (0) is xq = cxd. 
In a practical numerical calculation, which is always carried out in between finite ranges, Xq 
is assigned an arbitrarily large but finite value. But as Xq approaches infinity, both ^^^(O, Xq) 
and ipn{^,Xo) diverge, and a regularization process is needed to make sense of Eq. (j^). 

In the form of Eq. (|l|), taking into account of the advantage of its similarity in form to 
the first order equation, ipn{ci, x) can easily be solved, say, using the Dalgarno-Lewis method 
or logarithmic perturbation method ^^ to obtain, 

/■^ 1 /■?/ 

^n{a,x) = -%1)q{x) dy dz {a - Vn)^oiz), (8) 



^For purpose of illustration, we consider QM on a half line [0, cxo). 



where a and h are appropriate constants to satisfy the boundary conditions, in agreement 
with the resuhs of GR. 

From Eq. (|^), one recovers the universal functions F{x) (given as Eq. (15) in Ref. |Q) 



F{x) = -Mx)l dy^— rdzV{z)tljl{z). (9) 

Jb ip^{y) J a 

Together with Eq. (|]), one sees that the n-th order perturbation energy En is given by 

where the functional J{V, x) is given by 

J{V,x) ^ rdy^ r dzV{z)ijl{z). (11) 

Jo ^p^{y) Jo 

and the boundary condition at the endpoints has been taken care of appropriately. 

Next, we would like to point out that at least in the example of the ground state of the 
x^ anharmonic oscillator, the expansion of En in Eq. (|I0[) can be ill-defined because both the 
numerator and the denominator diverge as xq ^ c>o. This can be easily seen by combining 
the well-known results that Vn is of polynomial form in the Bender-Wu [^ x^ anharmonic 
oscillator and the mean value theorem. 

From the form of Eq. (0), one does not expect a priori that in the limit Xq — > oo, 
the ratio J{Vn, xo)/J{l, xq) becomes finite and xq independent even though numerically this 
comes out to be so. Hence to make sense out of Eqs. ([T0| ) and (0), a regularization procedure 
is in order. One can justify the numerical result obtained by assigning an arbitrarily large 
but finite value to Xq only after the result is regularized and the limit is proven to exist. 

The regularization procedure being proposed here is similar to the one we previously 
used in the extension of logarithmic perturbation theory to excited bound states in one 
dimension by appropriately mixing in the ghost state [0. For the zeroth order solution 
(unperturbed state), instead of using the square integrable eigenstate wave function ipo, we 
can mix in the non-square-integrable ghost state Xo by defining 

^o(a;) = ^o{x) + iaxoix), (12a) 



p{x) ^ ^l{x). (12b) 

and 

J,\S] = r dy^ r dzp{z)S{z). (12c) 

Jo p[y) Jo 



Note that in Eq. ( |12b| ), p{x) is the ordinary square of \l/o(a;), not |\l/o(a;)p. Then Eq. (|T0|) 
can be rewritten on firm mathematical ground as 

(|TT|) Now, we can show that the hmit in Eq. ([T3|) is the well-defined. This follows from 

MS] = - dyMy)My)S{y) 

a Jo 

/ dy^l{y)S{y) + ---, (14) 

Jo 



a Jo 



where ■ • ■ is a cr independent term. Upon substituting Eq. (|Tip into Eq. (|I^), we recover 

P I^dy^i{y)Vniy) P > , . u/ ^ w f\ fiK\ 

"^ Io-dyro{y) ^l dyMyWM^Pn^^iy) (i5) 

which is the ordinary QM result upon using Eq. (^) and (^). Hence, we have provided a 
rigorous justification of the SC method. It is interesting to note that we have also utihzed 
non-square-integrable functions through the ghost state mixing. 

To sum up, we see that the SC method correctly gives the perturbation energies, but as 
a ratio of two divergent quantities. We have regularized it through ghost state mixing and 
our final result is independent of the mixing parameter a. It is only after establishing the 
existence of the limit in Eq. (|T3p that we can accept the numerical convergence in Eq. (|TUp 
advocated in the SC method. 
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